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ABSTRACT 


The  Incompressible  flow  Inside  a  porous  cylindrical 
tube  having  steady  radial  Inflow  and  rotating  about  Its 
axis  In  a  sinusoidal  manner  is  obtained  through  solution 
of  the  complete  Navler>Stokes  equations.  The  nature  of  the 
tangential  velocity  profiles  that  result  and  the  magnitudes 
of  the  mean  pressure  drop  and  the  amplitude  of  the  pressure 
fluctuations  that  are  Induced  by  the  tangential  velocities 
are  discussed  in  some  detail. 
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1.  Introduction 

In  reference  1  (hereafter  referred  to  as  Part  I) 
we  considered  the  behavior  of  a  class  of  three-dimensional 
steady  solutions  of  the  Navler-Stokes  equations.  The 
flows  considered  were  those  such  that,  In  a  cylindrical 
coordinate  system  (r,  <t>,  z),  the  radial,  tangential,  and 
axial  velocity  distributions  were  of  the  form 

u  ■  u(r) 
v  »  v(r) 
w  »  zw(r) 

while  the  pressure  distribution  was  of  the  form 
p  -  Cjz2  +  px(r)  +  p2(r) 

In  the  above  expression  for  the  pressure  one  may 

p 

Identify  the  term  c^z  with  the  pressure  necessary  to 
sustain  the  axial  velocities,  the  term  p^(r)  with  the 
pressures  necessary  to  sustain  the  radial  motion,  and  the 
term  p2(r)  with  the  pressures  necessary  to  sustain  the 
tangential  velocities. 

This  special  class  of  steady  solutions  of  the  Navler- 
Stokes  equations  is  such  that  the  equations  for  the  radial 
and  axial  components  of  velocity  are  decoupled  from  the 
equations  governing  the  tangential  velocity  distribution 
and  the  distribution  of  pressure.  Because  of  this  decoupling 
of  the  equations,  one  finds  that  for  very  little  additional 
labor  one  can  obtain  the  behavior  of  a  simple  class  of 
unsteady  solutions  of  the  Navler-Stokes  equations ;  namely, 
solutions  of  the  form 
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u  *  u(r) 
v  -  v(r,  t) 
w  =  zw(r) 

p  =  c1z2  +  p1(r)  +  p3(r,  t) 

We  see  from  the  above  expression  for  the  pressure  that 
the  pressures  necessary  to  sustain  the  axial  and  radial 
velocities  are  identical  to  those  of  the  steady  flows,  as 
might  be  expected,  and  only  the  pressure  necessary  to 
sustain  the  tangential  velocities  p3(r,  t)  becomes 
dependent  on  both  position  and  time. 

As  in  Part  I  of  this  report,  these  flows  may  be 
visualized  as  the  flow  inside  a  cylindrical  porous  tube 
which  is  subject  to  a  rotation  about  its  axis  which  is 
some  arbitrary  function  of  time. 

In  what  follows,  we  will  report  in  some  detail  the 
behavior  of  these  solutions  when  the  porous  tube  is 
subjected  to  a  sinusoidal  oscillation  about  its  axis. 
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2.  Basic  Equations 

The  basic  equations  governing  the  fluid  motions  that 
we  wish  to  discuss  (see  Part  I,  Section  2)  are 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


i  -~(ru)  +  w  -  0 


r  dr 


U-—  -Iv2 

dr  r  p  dr  arlr  dr  J 

|»  +  »(„)] 

St  r  or  drLr  dr  J 


dw  -2 
zu~  +  zvr 
dr 


1  dp 

"p  fz 


1  d  f  dw  \ 
+  vz  i  r~) 
r  dr\  dr/ 


We  note  that  Equations  (2.2)  and  (2.4)  which  govern  the 
radial  and  axial  velocities  and  the  boundary  conditions  to 
be  applied  to  these  equations,  namely  u(0)  »  0,  u(R)  -  U, 
w* (0)  ■  0,  and  w(R)  *  0,  are  identical  to  those  discussed 
in  detail  in  Part  I.  The  radial  and  axial  components  of  the 
flow  are  therefore  identical  to  those  of  the  steady  flows 
of  Part  I. 

We  now  introduce  the  non-dimensional  variables 


(2.5) 

x  -  (r/R)2 

(2.6) 

t  =  St/R 

(2.7) 

f (x )  -  -2ru/RS 

(2.8) 

h(x,  t)  «  rv/RV 

where  S  and  V 

are  velocities  and  R 

a  length  that  will 

be  defined  later. 
(2.3)  becomes 

In  terms  of  these  new 

variables.  Equation 

(2.9) 

xh„„  +  nfh„  ■  nh_ 

XX  X  t 
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where  the  subscripts  Indicate  differentiation  with  respect 
to  the  variables  x  and  t  and  the  parameter  n  -  RS/4v 
is  a  Reynolds  number. 

Now  let  R  be  the  internal  radius  of  a  porous  cylinder 
oscillating  in  such  a  manner  that 

(2.10)  v(R,  t)  =  V  cos  cot 

Applying  this  boundary  condition  to  h(x,  t)  gives  us 

(2.11)  h(l,  t)  =  cos  At 

where  A  «  Rev'S.  We  also  require  that  the  tangential 
velocity  vanish  at  r  »  0,  so  that  we  have 

(2.12)  h(0,  t)  «  0 
If  we  write 

(2.13)  h  *  R. P.-^A(x)e”~^t|-  *  X  cos  At  +  u  sin  At 
where  A(x)  »  X(x)  +  ip.(x),  then  Equation  (2.9)  becomes 

(2.14)  xA"  +  nfA'  -  -inAA 


The  boundary  conditions  (2.11)  and  (2.12)  are 


(2.15)  A(0)  =  0  A(l)  -  1 


Let  0(x)  »  a(x)  +  i{5(x)  be  a  solution  of  Equation  (2.14) 
when  0(0)  «  0  and  0'(O)  «  1.  Then  the  solution  we 
require  is 


A(x) 


0(1) 
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or  writing  $(l) 
(2.16) 

(2.17) 

where  a(x)  and 


+  10^  we  have 


\(x) 


+  ^0 

a?  +  p2 


M-(x) 


axP  -  pjO 

a2  +  B2 
1  +  pl 


p(x)  are  solutions  of 


(2.18) 


xa"  +  nfa'  -  nflfJ 


subject  to  a(0)  =  0  and  a!(0)  =  1  and 
(2.19)  xp"  +  nfP*  *  -nfia 

subject  to  P(0)  =  0  and  ^,(0)  =  0. 

In  order  to  use  the  results  presented  in  this  report 
it  is  necessary  to  specify  the  characteristic  velocity  S. 
As  is  discussed  in  Part  I  (p.  38),  the  characteristic 
velocity  S  is  the  numerical  value  of  the  axial  velocity 
at  r  «  0  and  z  =  R.  The  exact  value  of  S  in  relation 
to  the  radial  inflow  or  outflow  at  r  *  R  is  obtained  from 
the  steady  solution  for  the  radial  and  axial  velocities. 
Figure  2.1  is  a  plot  of  the  ratio  of  n/|UR/v|  *  S/4|u| 
as  a  function  of  |UR/v|  for  the  simple  families  of  one- 
celled  flows  obtained  from  the  studies  reported  in  Part  I. 
As  is  indicated  in  this  figure,  one -celled  flows  under 
conditions  of  outward  radial  flow  exist  for  only  a  limited 
range  of  radial  Reynolds  number  (for  a  more  extended  dis¬ 
cussion  see  Part  I).  For  this  reason,  we  shall  limit  the 
computations  and  discussion  to  the  case  of  radial  inflows. 


'll 


Fig.  2.1  Relationship  between  the  characteristic  velocity 
S  and  the  numerical  value  of  the  radial  velocity 
U  into  a  porous  cylinder  of  radius  R  as  a  function 
of  the  Reynolds  number  UR/y  . 


6 


3 .  Velocity  Distributions 

An  analog  computing  machine  was  used  to  solve 
Equations  (2.18)  and  (2.19)  and  thus  determine  the  functions 
X(x)  and  u(x)  which  represent  the  tangential  velocity 
distribution,  for  a  range  of  both  the  frequency  parameter 
ft  and  the  Reynolds  number  n.  In  making  these  computations 
the  radial  velocity  function  f(x)  was  generated  by 
simultaneous  solution  of  the  appropriate  equation  for  f 
subject  to  the  proper  boundary  conditions  (Equations 
(2.2.6),  (2.2.8),  (2.2.9),  and  (2.2.10)  of  Part  I).  Typical 
tangential  velocity  profiles  at  an  instant  in  time  and 
envelopes  of  the  velocity  profiles  for  all  times  are  shown 
in  Figures  3*1»  3.2,  3-3»  and  3.^  with  the  Reynolds  number 
n  as  a  parameter  for  four  values  of  the  frequency  parameter 
(ft  «  0,  0.3,  1.0,  and  3.0  respectively).  In  Figure  3-5 
the  envelopes  of  the  velocity  profiles  are  shown  for  the 
particular  Reynolds  number  n  »  10  with  ft  as  a  parameter. 


at  o.*  0.6  ad  f.o 

r/R 


Fig.  3-1  Tangential  velocity  profiles  /I  =  0. 


Fig,  3.<J  Tangential  velocity  profiles  A  .  0,3, 


Fig,  3.4  Tangential  velocity  profiles  J\  *  3,0 


7 


4.  pressure  Distributions 

As  pointed  out  In  the  introduction,  the  preeaure 
distributions  for  the  flows  under  discussion  are  of  the 
form 

(4.1)  p  -  c^2  +  p2(r)  +  p3(r,  t) 

we  may  show  this  in  the  following  manner.  Prom  Equation 

(2.2)  we  note,  since  u  and  v  are  functions  of  r 
alone,  that  d2p/drdz  •  0.  Differentiating  Equation  (2.4) 
with  respect  to  r  gives 

C-21  ?  •  0  -  -*£[“§  +  7 

Prom  this  equation  we  see  that 

(4.3)  ^  +  w2 . 1  ^(r^)  -  constant 
or 

(4.4)  ||  .  Cl* 
and 

(4.5)  P  -  CjE2  p(r,  t) 

Since  in  what  follows  we  are  primarily  Interested  in 
the  effects  of  the  tangential  motions  on  the  pressure  we 
shall,  in  view  of  Equation  (4.5),  define 

(4.6)  .(X,  T)  .  P(Mi*) 


Zn  addition  we  define 
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(4.7) 

(4.8) 

(4.9) 

(4.10) 

(4.11) 


rU(jt)  -/ 


igdX 


▼lg(3t) 


-/> 


’,22(x)  */  ^ 

1  * 

»3  (x)  .  |[f*(l) 

T  4  (*)  *  *  ^jjf '  (*) 


Substitution  of  Equations  (4.6)  through  (4.11)  into 
Equation  (2.2)  yields 

▼  »  ^V2(T1icos2Qr+2T12sln  Or  cos  Ot+irggSln^x)  +  S^v^+v^) 

»  ^[(•‘n  •  t22^co*  20t  +  ^12  8ln  2^] 


+  +  *22*  +  S2(ir3  +  t4) 


(4.12) 


In  this  equation,  the  term  S2(t^  +  v^)  represents  that 
part  of  the  pressure  required  to  sustain  the  radial  motion, 
the  term  S2*^  representing  the  pressure  due  to  lnviscld 
effects  and  the  term  S2v^  representing  the  pressure  due 
to  vlseous  effects.  The  remaining  terms  on  the  right-hand 
side  of  Equation  (4.12)  represent  the  pressure  required 
to  sustain  the  tangential  velocities.  There  are  two  effects* 
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a  mean  pressure  difference 

(4.13)  —5  -  |(*n  +  *22* 

and  a  fluctuation  in  pressure  at  twice  the  frequency  of  the 
oscillation  of  the  porous  boundary  which  has  an  amplitude 


(4.14) 


wns. 


-  T, 


22 
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Figures  4.1  and  4.2  are  plots  of  these  two  quantities 
evaluated  at  r  •  %  »  0  as  obtained  from  the  analog 
solutions  for  X  and  u  discussed  in  Section  3  for  the 
Reynolds  number  range  0  <  n  £  55  end  for  values  of  the 
frequency  parameter  equal  to  0,  0.3,  1.0,  and  3.0.  For 
purposes  of  comparison,  there  is  also  plotted  in  Figures 
4.1  and  4.2,  the  mean  pressure  drop  and  amplitude  of  the 
pressure  fluctuation  (which  in  the  limiting  case  fl  0 
is  equal  to  the  mean  pressure  drop)  obtained  from  a  form  of 
Burgers'  analytical  solution  (Reference  2)  for  the  steady 
(fi  a  0)  tangential  motion  which  is  valid  for  large  n, 
namely 


Apm  and  APf 


(4.15) 


v 

V 


1  1  -  e 


-nx 


Vx  1  -  e‘n 

If  Equation  (4.15)  is  used  to  obtain 
for  fl  0,  there  results 


ipv2 


l**L» 

ipv5 


|in»  -i 
2  2 


(n»l) 


S  /.#  s  *  /o  to  SO  /oo 

r> 


Pig.  4.1  Mean  pressure  drop  in  an  oscillating  vortex  tube 
as  a  function  of  Reynolds  number  for  various  non- 
dimensional  oscillatory  frequencies  fL  of  the 
vortex  tube . 


.■I.— . L .  ^ _ _ J _ L. . . . . . -J. 

•S  10  c  s  to  to  so 

n 


Pig.  4.2  Amplitude  of  the  pressure  fluctuation  in  an 

oscillating  vortex  tube  as  a  function  of  Reynolds 
number  for  various  non-dimensional  frequencies  A 
of  the  vortex  tube.  The  non-dimensional  frequency 
of  the  pressure  fluctuation  is  2 A  . 
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5.  Discussion  of  Results 

It  Is  convenient  to  start  the  discussion  of  the 
results  Just  presented  with  an  analysis  of  Plgures  4.1  and 
4.2.  We  note  that  for  n  -b  0,  4pm(0)/^pV2  »  -0.5  and 

|Apf (0)  1,^,,,/ipV2  *  0.5.  This  is  easily  shown  to  be  true 
because,  in  the  limit  of  zero  Reynolds  number  or  infinite 
viscosity,  the  fluid  within  the  vortex  tube  turns  as  a 
solid.  Thus  as  n  -b  0 


(5.1) 


v  -  v(R)g 


and,  from  Equation  (2.2),  that  portion  of  p  that  is  due 
to  the  tangential  velocity  is  then  obtained  from  an  integra¬ 
tion  of 


(5.2) 


*E.£XI 

dr  r 


Substituting  Equation  (5.1)  into  (5.2)  and  integrating 
yields 

(5.3)  P(0)  -  p(R)  --4p[v(R)]2 

Since  we  have  chosen  v(R)  »  V  cos  t  we  have 


(5.4)  p(0)  -  p(R)  «  -ipV2  cos2<nfc 

or 

■  -  4(1  +  COS  2oot ) 

Thus  for  n  -b  0  the  mean  part  of  the  non-dimensional 
pressure  difference  is  -0.5  and  the  amplitude  of  the 
non-dimensional  pressure  fluctuation  is  0.5.  The 
frequency  of  the  pressure  fluctuation  is  twice  the  frequency 
of  oscillation  of  the  boundary. 

Whether  or  not  the  mean  pressure  drop  and  the  amplitude 
of  the  pressure  fluctuations  become  larger  numerically  as 
the  Reynolds  number  n  is  Increased  depends  on  the  frequency 
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parameter  0.  For  frequency  parameters  smaller  than  unity, 
the  general  trend  le  for  the  pressure  difference  and  ampli¬ 
tude  to  Increase  numerically  as  n  increases.  As  0 
becomes  larger  than  unity,  the  tendency  is  for  the  mean 
pressure  differences  and  amplitudes  to  at  first  decrease 
numerically  as  n  becomes  larger  than  unity  and  then 
finally  Increase  at  some  Reynolds  number.  The  Reynolds 
number  at  which  the  pressure  fluctuations  start  to  Increase 
is  a  function  of  ft.  The  larger  ft  becomes,  the  larger  is 
the  range  of  n  above  unity  for  which  the  mean  pressure 
differences  and  the  amplitudes  of  the  pressure  fluctuations 
due  to  the  oscillating  boundary  are  of  small  magnitude. 
Eventually,  for  an  incompressible  fluid,  all  curves  such 
as  those  shown  in  Figures  4.1  and  4.2  must  turn  upwards 
because,  given  any  finite  frequency  no  matter  how  high,  one 
can  always  choose  a  Reynolds  number  high  enough  so  that  the 
shears  developed  between  successive  maxima  of  the  velocity 
distribution  curves  become  negligible.  When  this  condition 
exists,  the  fluid  Introduced  through  the  porous  boundary 
is  able  to  maintain  to  some  degree  its  initial  angular 
momentum  with  the  result  that  high  tangential  velocities  are 
carried  toward  the  axis.  In  an  incompressible  medium  these 
high  tangential  velocities  cause  pressure  fluctuations 
that  are  transmitted  Instantaneously  to  the  boundaries  of 
the  flow.  It  appears  from  what  has  been  said  above  and 
from  the  behavior  of  the  solutions  which  have  been  obtained, 
that  for  all  finite  frequencies  the  mean  pressure  difference 
and  the  amplitude  of  the  pressure  fluctuations  become 
proportional  to  n  as  n  ->  «,  as  has  been  discussed  for 
the  case  ft  ■  0  in  Section  4. 

The  reason  for  the  initial  drop  in  the  numerical 
values  of  Apm(0)  and  |Apf (0) |  as  n  Increases  for 
frequency  parameters  greater  the  unity  can  best  be  explained 
as  follows.  At  low  Reynolds  numbers  (n  «  1),  the  viscosity 
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of  the  fluid  le  »o  high  that  reasonably  large  tangential 
velocities  are  maintained  all  the  way  to  the  center  of  the 
vortex  tube  by  strong  vlsoous  action.  As  the  viscosity  is 
decreased  the  fluid  toward  the  center  of  the  tube  no  longer 
feels  so  strongly  the  effect  of  the  surrounding  cylinder  at 
r  »  R.  However,  the  vlseoslty  Is  suoh  that  the  entering 
fluid  can  not  retain  its  Initial  angular  velocity.  The 
result  Is  that  the  velocities  near  the  oenter  of  the  tube 
are  not  as  large  as  they  are  for  n  -  0  as  may  be  seen 
from  the  tangential  velocity  envelope  for  0  •  3.0  and 
n  -  10  shown  In  Figure  3.5.  It  should  be  obvious,  for  a 
fixed  Reynolds  number  n,  that  the  velocity  profiles  and 
their  envelopes  should  become  somewhat  similar  to  those  In 
the  vicinity  of  an  oscillating  flat  plate  (see  Reference  3) 
as  0  Is  Increased  to  larger  and  larger  values. 

The  reason  for  the  slight  discrepancy  between  the 
behavior  of  the  mean  pressure  drop  and  the  amplitude  of  the 
pressure  fluctuations  obtained  from  the  analog  computations 
and  the  analytical  results  obtained  from  Burgers'  solution 
for  0  »  0  was  a  first  thought  to  be  computational  Inaccuracy. 
After  this  possibility  had  been  checked  and  the  computational 
Inaccuracies  found  to  be  far  less  than  the  amount  of  the 
discrepancy.  It  was  realized  that  the  discrepancy  Is  due  to 
the  different  radial  velocity  profiles  for  the  two  cases. 

In  the  case  of  Burger's  solution,  the  radial  velocity  Is  of 
the  form  u/u(R)  -  V5F,  while  for  the  cases  Investigated  In 
this  report,  the  appropriate  limiting  form  of  the  radial 
velocity  distribution  as  n  •  Is 


=  ,ln  <T> 


(see  Part  I,  page  64).  It  Is  this  difference  in  radial 
velocity  distributions  that  is  responsible  for  the  differen¬ 
ces  in  the  two  curves  for  0  ■  0  and  n  •  in  Figures  4.1 
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and  4.2. 

It  Is  lnterastlng  to  note  that  for  the  particular 
caae  of  flows  for  which  the  radial  velocity  at  r  »  R  Is 
proportional  to  the  tangential  velocity  at  r  »  R,  the 
pressure  fluctuations  In  the  flow  will  vary  as  for 
high  enough  Reynolds  number  since.  In  this  case,  the 
Reynolds  number  n  Is  proportional  to  V  and,  for  high 
enough  Reynolds  number  for  a  fixed  frequency, 

|Apf (0)  becomes  proportional  to  n. 
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